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'. ABSTRACT. 
p ,. 

. In the early eighties, Lopez, Gonzalez-Arroyo and the present author 

^< ' proved that, if at a given Qq large enough for perturbative QCD to be valid, 

■ structure functions behave as a power of a; for x ^ 0, then for all larger 

one has 

g : FGix,Q^) ~ BG[asiQ^)]-^^x'^ 

^' Rix,Q^) = '-^^^, 

TT 

with Dii, d+, Bg, To calculable in terms of B^, A. Moreover, it was suggested 
that the "hard" part of the scattering cross section for real photons (Compton 
scattering) obeys a similar law, so that 



with a value of A comparable to that in the expression for the structure func- 
, tions, and where ~ log^ s is a universal, Pomeron-type cross section, and 

^ • Ajp, Bryp are constants. In the present paper it is shown that the recent HERA 

measurements may be described by these formulas, with a chi-sqared/d.o.f. 
substantially less than unity, and with values of the parameters compatible 
with those of the old fits of the '80s. Moreover, further discussions are pre- 
sented both on the low limit, and the transition between Compton and 
deep inelastic scattering, in particular in connection with possible saturation 
of the coupling constant as{Q^) at small and on the ultra high energy 
limit, and how one might test the so-called BFKL conjecture. 



lim F2{x,Q'^)^x' 



With respect to the last we find some evidence against it, at least at the HERA 
energies. 
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1 -INTRODUCTION 



We consider in this note high energy inclusive scattering of (virtual or real) photons off protons: 

7*('?)+p(p)^all. (1.1) 

In the case of real photons we will call the process Compton scattering; for virtual photons we have DIS 
(=deep inelastic scattering). In this last situation we consider the region of very small x with* 

x = Q^/y,y = p-{p + q),Q'' = -q^. (1.2) 

Compton scattering may be thought of as the limit of DIS when — > 0, a; ^ in such a way that 

Q'^/x^s, s = {p + qf. (1.3) 

To a large extent the present work may be considered as an aggiornamento of the old analysis of C. Lopez, 
A. Gonzalez- Arroyo and the present author, applying it to the recent HERA DIS and Compton data. In this 
respect it is shown that the new data are fitted very well by the formulas derived from QCD in the small x 
region, and that the expressions 

F2{x, Q2) ~ B5[a,(Q2)]-d+^-A + B^^[a,(Q2)]-D.,^o.5^ 

Fg{x,Q^) ~ BG[as{Q^)]-'^x-^ (1.4) 
R{x,Q-) = '-^^, 

TT 

with Bg, -Dii, d+, ro calculable in terms of Bs, A, give an excellent approximation to the data for x ^ 10~^ 
(Sect. 2). Likewise (Sect. 3) we find that the very high energy Compton cross section is still correctly described 
by the formula 

c^p — B^pS -\- A^pCr , (1.5) 

with a value of A similar to the A in (1.4), and where ~ log^ s is a universal, Pomeron-type cross section. 
Ayp and B^p arc constants. 

The quality of the fits is so good, with chi-squared of less than one by d.o.f., that one may consider 
studying the corrections to (1.4), and the interpolation between (1.4) and (1.5); this we do in Sect. 4. In 
Sect. 5 we discuss the ultra high momentum limit of our formulas, in particular in connection with the BFKL 
approach which suggests the asymptotic behaviour 

F2{x,Q^)r~.x'='"'\ (1.6) 

We study the region of validity of (1.4), limited for fixed a; by a certain Q'^{x). We get indications that, for 
x < 10-2, Q^{x) - 200 - 300 GcV^ 

In what regards the transition from (1.4) to (1.6) we get (not conclusive, however) evidence against it 
at least in the region covered by the HERA data, i.e., below ~ 800 GeV^. 



More details on notation may be found in ref.l 
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2 -DEEP INELASTIC SCATTERING AT a; ^ 
2.1. -General considerations 

For DIS the relevant quantities are the structure functions, 

F,{x, Q% Fg{x, R{x, Q') = F^i^^Qp - Fdx,Q') _ (2.I) 

Fq is the gluon structure function, and R is (proportional to) the longitudinal one. 

It is convenient to use the singlet function Fg, normalized so that the momentum sum rule reads 



i 



1 

2\ , 77 („ r>2^ 



dx[Fs{x,Q^)+FG{x,Q^)] 1. (2.2) 

The relation between F-2, and Fg is as follows: one has 

F2(x, = (e2) \Fs{x, Q^) + Fj,s{x, Q^)] , (2.3a) 

where Fj^s is the so-called nonsinglet structure function. Since in the limit a; ^ it decreases very fast 
compared to Fs, we may consider the approximate relation 

F2(x,Q2) ~ (e2)F5(x,Q2). (2.3b) 

x—>0 ^ 

The average charge of the excited flavours in (2.3) is (e^) = ^ for n/ = 4 or (Cg) = for n/ = 5. For the 
values of 0^ a; with 8 GcV^ < < 65 GeV wc arc in a mixed situation in which bottom is excited in the 
s variable but not in the variable. We will thus present results for both Uf = A and 5: as will be seen 
they differ very little. For ^ 9 GeV^ only three flavours are excited in the Q"^ variable. Both for this and 
other reasons this energy region requires a specific treatment. 

The evolution equations for the structure functions may be written in two ways. One has the so called 
Altarelli-Parisi, or DGLAP equations ^'^1, that we write only for the singlet structure functions, 

dF,{x,Q^) _ ^ , , d^p..(^)^.(^/^^g.)^ (2.4) 



dt 



and the explicit form of the splitting functions Py may be found in refs. 1,3. Alternatively, one may define 
the moments of the structure functions, 

pis{n,Q^) = jldxx^-^Fs{x,Q^), 
^lG{n,Q^) = J^dxx^-^FG{x,Q^) 

and then the integro- differential equations (2.4) may be solved for the /Xj: 



D(«) 



(2.6a) 



where the anomalous dimension matrix D(n) is 



[4] 
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2n{n + 1) 

^2 I 



D(n) = 



+ 



16 

33 - 2n/ 

3n/ n2+n + 2 
~8~ n(n + l)(n + 2) 



V 2n(n2-l) 4n(n-l) ' 4(n+l)(n + 2) 

The function is related to the digamma function, 



+ 



33 -2n/ 951(71) 
16 



(2.6b) 
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Si{z) = i^iz + 1) + 7B, ib{z) = dlogr(z)/dz, 7B = 0.5772 .... 

It should perhaps be stressed that eqs. (2.4) and (2.6) are strictly equivalent, as known from the very first 
works on the subjectl^l; (2.6) follows from (2.4) by taking moments and integrating, and (2.4) from (2.6) by 
inverting the Mellin transform (2.5) and differentiating. 

Eqs (2.6), in the limit a; — > 0, were solved long ago by Lopez and the present author for Fjvs, F2 and 
Fg to leading order (LO) in ref.5; including next to leading corrections in ref.6 and, for the longitudinal 
structure function by the quoted authors and Gonzalez- Arroyo in ref.7 [where the interested reader may find 
the extension of eqs. (2.4), (2.6) to R{x, Q^)]. The comparison with experiment was carried to leading order 
in refs. 5,7 ; the comparison with experimental data of the calculation to next to leading order (NLO) was 
given in ref.8. 

The LO results are remarkably simple. It was proved in ref.5 that, under certain conditions, the structure 
functions are given by the very explicit expressions, which therefore solve the evolution equations, * 



Fs{x,Q^) Bs[a,{Q')] 

X— »0 



-d+(l+A)^-A 



Fg{x,Q^) 



a:->0 



Fs{x,Q^), 



x— >o 

R{x,Q^) ^ ro(l + vl)^^. 

X— >0 TT 



We will also use the quantities 



Bs 



{^q)Bs, BjMS 



ie^B^S- 



(2.8a) 

(2.8b) 
(2.8c) 
(2.8d) 

(2.8e) 



A, Bns and Bs are free parameters, independent of x, Q^, although they may have a dependence on n/ 
(expected to be slight; it will be discussed in Sect. 6). p is known from Regge theory to be the intercept of 
the — p trajectory, p ~ 0.5. The quantities B^ro are obtainable in terms of A, Bg: 



Bail + X) d+(l + X)-Dn{l + X), 



Bs 



ro(l + A) = 



3(2 + A) 



1 + 



^12(1 + A) 
3n/ d+(l + A)-i:'ii(l + A) 



(2.8f) 



2 (3 + A)I)i2(l + A) 
Here the Dii(l - p), Ai(l + A) and the d±{l + A), 

d±{l + A) = i {£)ii(l + A) + £)22(1 + A) ± ^/[Du{^ + A) - £>22(1 + X)? + 4£'i2(l + A)£'2i(l + A)} 

are the matrix elements and eigenvalues of the matrix D(n) evaluated (respectively) atn = l — p, n=l + A 
(the eigenvalues ordered so that d+ > d_). We note that eqs. (2.5), (2.6) are valid for arbitrary (even 
complex) A-alues of n. Thus. (^qs. (2.8) give an explicit expression for all three structure functions Fs, Fq 



* For ease of reference we will reproduce the proof in Sect. 6 here. 
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and R in terms of the only two free parameters A and Bs (-Bats also intervenes if we include Fns in the 
analysis) . 

Surprisingly enough, eqs. (2.8) and, more generally, the work of refs.5-8 seem to have been ignored by 
recent publicationst^'"'^''! . Even more surprising, the HERA physicists have not used eqs. (2.8) in spite of 
their simplicity to analyse the experimental data, relying instead on a painstaking numerical integration of 
the Altarelli-Parisi equations (2.4) -a procedure not only infinitely more complicated than use of the explicit 
solutions (2.8), but that, as will be shown, produces larger errors. 

2. 2. -Low and medium x analysis, and higher and very small x predictions 

In the 1980 paper, Lopez and the present author analyzed DIS data in the range 

3.25 GeV^ < < 22.5 GeV^ x>2x 10"^. (2.9) 

We found that the parameter A could be well determined while Bs (which gives the overall normalization) 
was more uncertain. This uncertainty was due basically to the strong contribution of subleading effects, 
especially those due to the nonsinglet structure function, at the (relatively) large values of x used. One had 

A = 0.37 ± 0.07, 0.02 > Bs > O.OOL (2.10a) 

[Actually the value and error take into account also the analysis of R; F2 alone would produce a slightly 
smaller A, A = 0.36 ]. The results of the evaluation were subsequently shown to be essentially unaltered 
by inclusion of NLO corrections!^! . As for B^s, the fits based on functions containing both singlet and 
nonsinglet parts fix it only within large errors. Fortunately, however, the W3 structure function in neutrino 
scattering is purely nonsinglet and allows a reasonably precise determination: B^s ~ 1.4 to 1.8 orl^l, 

Bjvs ^ 0.3 to 0.6. (2.10b) 

With these numbers it is a trivial matter to predict the very low x measurements of at HERA in terms 

of the only parameter Bs that has to be taken as essentially free because, as explained before, it is not well 
fixed by the old, larger x data. In a more precise determination we will of course also allow A, Bpfs and even 
A to vary. 

We will perform the analysis in two steps. In the first we will consider the old, 1993 datal^^l, whose 
errors were typically 0(10%), and with points in the range x < 10^^,8.5GeV^ < < 65GeV^; discussing 
then the more recent, as yet unpublished set of data which have much smaller errors, and cover a wider 
range. Later in Sects. 4, 5 we will consider the theoretical corrections to the formulas used. 

The results of the analysis of the 1993 data are presented in the following two tables. There the value 
of A choosen is 0.38, as it gives the best fit. The approximation of neglecting Fns, i.e., Eq. (2.3b) was used. 
For the QCD parameter A we first take it such that we reproduce the value of asint^) = 0.32, i.e., 

A{nf = 4, 1 loop) = 0.200 GeV, A{nf = 5, 1 loop) = 0.165 GeV. 



A d+(l + A) Bs Bg{1 + \)/Bs ro(l + A) xVd.o.f. 

0.38 ±0.01 2.406 ±0.100 (2.70 ± 0.22) x IQ-^ 20.56 ±0.54 6.24 ±0.24 9.13/(32-2) 

Table la.- n/ = 4; .4(1 loop, n/ = 4) = 0.200 GeV; as{ml) = 0.32. 



A d+(l + A) Bs Bg{1 + \)/Bs ro(l + A) xVd.o.f. 

0.38 ±0.01 2.595 ±0.110 (2.19 ± 0.11) x lO'^ 18.53 ±0.57 6.20 ±0.29 9.34/(32-2) 

Table 16.- n/ = 5; yl(l loop, n/ = 5) = 0.165 GeV; as{ml) = 0.32. 
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Fig.l.- Comparison of theoretical prediction, eq. (2.11) with experiment. 

One can consider fitting also the QCD parameter, A. In this case one discovers that, due to the slow, 
logarithmic variation of as with A and the large size of the experimental errors, the effect of altering A may 
be largely compensated by a change in Bs- As an indication, we give the results of an evaluation with a 
small A: 



A rf+(l + A) Bs Bg{1 + X)/Bs ro(l + A) xVd-o.f- 

0.35 ±0.01 2.735 1.19 x lO'^ 22.53 6.93 7.88/(32-2) 

Table II.- n/ = 4;yl(lloop,n/ = 4) = 0.10 GeV; as{ml) = 0.28. 

The difference between these results and those of Tables I a, & are to be considered as an added uncertainty 
in all that follows. For the sake of conciseness, we will however consider most of the time only the parameters 
reported in Table la. The graphic comparison with experiment is shown in Fig.l. Here, because of the 
smallness of the x^) "we have only given the curve for the central values, 

FM^) = (2.70 X 10-») |„.(g^)l-»x-«« .,,(0^) = (33_,„^)i„^'^'.^(„,,GeV)^ - <^-"" 

Morover, the error bars represent the statistical plus systematic errors composed quadratically. 

We find the phenomenon already encountered in the analysis of the old, low energy data: A is well 
determined, but Bs is less precisely fixed. Nevertheless, the errors have diminished substantially, and the 
consistency between the present analysis and that of refs. 5-8 is, to say the least, remarkable. 
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Figure 2.- Prediction for the gluon struc- 
ture function with n/ =4. For = 20 
GeV^ also the evaluation with n/ = 5 is 
shown. The points at that value of 
arc the results of the DGLAP calculation, 
with statistical errors only. 



The results obtained from the fit to F2 permit us to calculate also, and without any extra parameter, the 
two remaining structure functions, which for the central values of the parameters are, 



Fg(x,02) = (55.5x10-3) K(g2)] 



-2.406^-0.38 



R = 6.24 



(2.11b) 



(2.11c) 



with as{Q^) still as in (2.11a). 

We do not plot R as there are unfortunately no experimental data available. The gluon structure 
function is depicted in Fig. 2, where for the sake of comparison we also show the results based on a numerical 
evaluationt^^l using experimental data for F2 and the DGLAP equations with the method of ref.lO. The 
superior accuracy and simplicity of our formulas should be evident: note that the largest source of error in 
our evaluation is the uncertainty in the number of flavours between 4 and 5. 
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2.3.- Tie new HERA results. 



We will here consider the comparison of our predictions with the more recent HERA data (as yet 
unpublished). If one replaces, tels quels, the formulas used for the fits to the old data, then a very large 
chi-squared is produced. The reasons for this, rather obvious, are two. First of all, the precision of the data 
is such that the approximation of neglecting the F^s piece is no more justified. Secondly, the range is now 
such that the one loop evolution of is not sufficiently accurate. If we restrict our analysis to a; < 10~^, we 
have 63 data points with varying in the bounds 12 GeV^ < < 350 GcV^. We do not consider in the 
same fit smaller values of because this would imply a large variation of n/, on which A, Bs are expected 
to depend, and because the NLO corrections would certainly play an important role. 

We consider two possibilities: a restricted fit, only in the region 

12 GeV^ < Q2 < 90 GeV^ 

where it is sufficient to include the NS contribution, and an evaluation in the full range where the two loop 
expression for has to be taken into account as well. 



A Bns Bs ro x^/d-o.f. 

0.39 ±0.01 1.03 2.73 x 10-3 6.03 44.2/(48-3) 

Table III. - n/ = 4; a; < lO'^; 12 GeV^ < < 90 GeV^ ; A{1 loop, n/ = 4) = 0.2 GeV 



The results of the fits are again excellent. For the first case they are summarized in Table III, where 
we also give the parameters pertinent to the NS contribution, and A is considered fixed. For the full range, 
the results are given in Table IV. Note that the value we obtained of A, taken now as a free parameter in 
the fit, cannot be considered as a true determination, as we have not included second order effects (that we 
discuss below). 



A yl(21oops,n/ =4) d+(l + A) Bns Bs ro x^/d.o.f. 

0.38 ±0.01 0.10 ±0.01 GeV 2.406 0.772 (1.0 ± 0.22) x lO'^ 6.24 50.41/(63-4) 

Table IV. - n/ = 4; a; < lO'^; 12 GeV^ < < 350 GeV^ 



We do not give the ensuing expressions for Fg or R. If we take the values of the parameters, and the range 
of the variables pertaining to Table III, they are like those in eqs. (2.11), mutatis mutandi. If we consider 
the situation in which we take as to two loops. Table IV, then the expressions for Fq, R would not be more 
reliable than before because the two loop corrections to these quantities differ from the ones to F5. Thus 
we give only the expression for the structure function F2. We have, 

F2{x, Q2) = (1.0 ± 0.22) X 10-3[a,(Q2)]-2.406^-o.38 ^ o.77[a«(02)]-o-5i4^o.5^ ^2.12a) 

where now 

(n^\ = l^TT / _ 153 - 19n/ logloggyyl^ 
"^'^^ ' (33-2n/)\ (33-2n/)2 logQV^' 

The comparison with experiment of the results from the restricted fit using the values of the parameters 
recorded in Table III woud give something very much like one what sees in the Figures 1 and 2. The 
comparison of the theoretical evaluation for F2 using eq. (2.12) and experiment is shown in Fig. 3. 



j, n/ = 4, yl = 0.10±0.01. (2.12b) 
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Figure 3.- Fit to the latest 
HERA data with eq. (2.12) 
and Q2 from 12 to 350 GeV^. 
The errors of the experimental 

points, of the order of the size 
of the dots, are not shown. 



To end this section we fit separately the low data points, with a two loop expression for . We then 
find the results of Table V. 

A yl(21oops,n/ = 3) B^s Bs xVd-o.f. 

0.29 ±0.02 0.14 GeV 0.75 0.3 x lO'^ 18/(45-3) 



Table V. - n/ = 3; x < IQ-^; 1.5 GeV^ < < 8.5 GeV^ 



We do only give errors on A. In fact the results in Table V should be taken as representing effective 
estimates. The reason is that the values of are so low that an exact treatment of subleading corrections 
would be essential to get realistic values of the parameters. This includes not only 0{as) corrections, but 
diffractive ones as well. This may easily alter substantially the results reported here; one should realize, for 
example, that even without the corrections a chi- squared/ d.o.f. of less than one may be obtained provided 
A ^ 0.35, with Bs ^ 0.6 x 10~^ and .Bjvs ~ 1-1. Note also that the compatibility of the results reported 
in Tables IV and V is made more apparent if we compare values of the Bs, B^Si i-C-) we extract the factor 
(el) = 5/18 (n/ = 4) or 2/9 ( for n/ = 3). Then we get, 

Bs = 3.6 X 10"^ (Q^ > lOGeV^); (1.4 to 2.7) x 10-^(Q^ < lOGeV^); 
Bns = 3.11 iQ^ > 10 GeV^); 3.4 to 4.4 (Q^ < 10 GeV^). 
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3 -HIGH ENERGY COMPTON SCATTERING 



The results of the previous section indicate that a behaviour 

F2{x,Q^) ^Js{Q^)s^ + fNsms'' + ... (3.1) 

x—>0 

where s is the cm. energy squared, and we set the scale so that it is measured in GeV^, may hold beyond 
the region of applicability of perturbative QCD; note that s ~ Q^/x. This was first suggested in ref.l2 for 
Compton scattering, 7 + p — > all. The supposed subleading terms, denoted by dots in (3.1) were identified 
with diffractive effects the Pomcron. To be precise, one may suppose that the photon has a "hard" , or 
pointlike component to which (3.1) applies (without the dots); and a "soft", hadronic or Pomeron component, 
which in old fashioned photon physics'^^l was identified as connected with the probability of finding the rho 
resonance in the photon. One thus writes 

a^p{s) = C7^j,(s) + CT^p(s), (3.2a) 

where 

a^p(s) = B^p{s/1 GeV^)^, (3.2b) 

a^p{s) = A^pa^{s) (3.2c) 

and ^^p{s) is a universal, Pomeron hadronic cross section [into which, for reasons of convenience, we have 
also added the p — /° trajectory contribution, corresponding to the term fNs{Q^)s^'^ in (3.1)]. 

In view of its universality, <5"^(s) may be obtained, up to a constant, from any hadronic scattering cross 
sections, say 

C7rp(s) = CttpCT (s); (Tp_p;j,(s) = Cp—p-pa (s). i^-^) 
The formulas employed for fitting the various hadronic cross sections are thus, 

o-hpis) = Chpa{s), (3.3) 

where the constant Chp depends on the process and a{s) is obtained saturating the Froissart bound in the 
improved version of ref.l4* (but with the overall constant and an additive constant left as free parameters), 
plus a Regge pole contribution corresponding to the /° trajectory, which as explained we have found it 
convenient to incorporate into a^: 

a{s) = A^log-—4—— + l+Afo{s/lGeVY'-'- (3-4) 

mi log 2 (s/m4) 

The values of Ap, Ap and Chp are obtained fitting the cross sections ctttp = ^[o'-K+p + fir-p] ^-nd a^p = 
\Wpp + <^pp]- For energies above 500 GeV we assume app = app + "Regge", where the piece "Regge" 
is obtained extrapolating the low energy difference app — app to high energy with a Regge pole formula, 
o'pp — CTpp ~ Cs^p. This gives a minute correction, 0.6 to 0.2 mb, for ~ 500 to 2000 GeV. One then finds 
the universal parameters 

Af = 0.0116, Afo = 0.69 (3.5a) 

and moreover 

C^p = 23.0 mb; Cjvp = 39.6 mb. (3.5b) 

Note that the ratio CNp/C-^p = 1.72 is reasonably close to the value 3/2 predicted by the naive quark model. 

It should be stressed that the present paper is not about hadronic cross sections. If a set of experimental 
measurements of 'JTp,pp,pp existed from "low" energies (s^ ~ 10 GeV) to the large values reached in 'yp 
scattering, 200 GeV, a fit like that in eqs. (3.4,5) would be unnecessary; and indeed, one could also 



* This corrects an error in ref.l4, where the scale factor is wrongly given as log'^(s/m,r) instead of the 
correct value log' (s/m^), cf. (3.4) below. 
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have used the phenomenological interpolations provided 
because of a theoretical prejudice in its favour, which is 
two parameters, shown in Fig. 4. 



by the Particle Data Groupli'"^!. Wo choose (3.4) 
substantiated by the quality of the fits, with only 
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Figure 4.- Fit to the up (indented) 
and WPp+ pp] cross sections with 
eqs. (3.4) and (3.5). 



With this we write a^p as in eqs. (3.2). The values of the constants one obtains are, for a x^/d.o.f. 
11.3/12 (with statistical errors only). 



C-yp = 102.0 ± 2.0 jjb; B^p = 1.64 ± 0.30 /xb; A = 0.23 ± 0.03. 



(3.6) 



The ratio C^pjC^p = 225 is close to the value 220 obtained at low energy from experiment and vector 
meson dominance arguments'^^l. Finally, the parameters A = 0.23 and B-yp = 1.64 are reasonably similar 
to those obtained at "low energy" (s2 < 20 GeV) in ref.l2, that is, A = 0.40 and B^p = 1.40, while A is 
also close to the value deduced in the previous subsection (2.3) from small x in deep inelastic scattering, 
A ~ 0.3. One should not take the small discrepancies very seriously; the values recorded above for the 
parameters were obtained taking only statistical errors into account in the recent, very high energy HERA 
datalis]. If we include also systematic errors, then we get A = 0.26 and a chi-squared/d.o.f. less than one 
provided 0.22 < A < 0.29, perfectly compatible with what we found before, modulo the (expected; see Sect. 6) 
dependence of A on the number of excited flavours. 
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Figure 5.- Fit to the 7p 
cross section with eqs.(3.2) 
In the indent, a blowup of 
the low energy region. 



The quahty of the fit should be apparent from Fig. 5, where only statistical errors are shown for the 
high energy HERA data; also impressive is the stability of the extrapolation by one order of magnitude in 
energy from the 1980 analysis to the present one. 
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4 -SUBLEADING CORRECTIONS 



4.1. -Two loop QCD corrections 

The equations we have used take only into account LO (one loop) QCD, except for the running of 
as{Q^). A proper analysis should fully take into account NLO (two loop) corrections. These are known 
since the work of ref.6, and amount to replacing (2.8a-d) by 



Fs{x,Q^) ^ Bs 

X— >0 



1 + 



vsil + X)asiQ^ 



47r 



-d+{l + \) 



X 



(4.1a) 



Bg{1 + \) 



1 + 



VG{l + \)as{Q^) 



47r 



Fs{x,Q^), 



Fns{x,Q^) c^^Bns 



1 + 



47r 



(4.1b) 



(4.1c) 



where the Vi (n) are known, complicated combinations of one- and the two-loop anomalous dimension matrix 
elements and the one loop Wilson coefficients whose explicit expression may be found in ref.8. Because 
of this structure, the bulk of the NLO effect may be approximated for a single structure function, by just 
leaving yl as a free parameter -precisely as we did in subsect. 2.3. The full NLO analysis, with eqs. (4.1) 
taken into account, requires more precision, in particular incorporating also other subleading effects and will 
be presented in a separate paper. 

4. 2.- The joining of low and high 

In Subsect. 2.3 and Sect. 3 we have shown that on can reproduce the cross sections 



(77*(-Q2)p(s) 

for, respectively, off shell (but with fairly small Q^, down to 1.5 GeV^) or on shell photons, = 0, with 
formidas with similar functional form for the structure function and the "hard" piece of the Compton cross 
section. It is then tempting to try and join the two approaches. This presents a number of problems, 
and rises some very interesting questions, that we now briefly discuss, leaving the full analysis for a future 
publication. 

First of all, it is clear that an equation like (2.8a) (say) cannot hold for too small. Not only, if taken 
literally, it would become singular for ^ A^, but, from the relation between F2 and (7^.(_Q2)p(s) it follows 
that, as Q'^ 0, 

'^.f^F,{x,Q% (4.2) 



(T, 



7*(Q=)P 



so F2{x,Q'^) must vanish proportional to for small. 

There are three ways in which one may look at eqs. (2.8), (4.2). First of all comes the question of the 
joinig of the two. If we take literally the values of the parameters we found in Sect. 2, and the hard part 



of the photon-proton cross section, this occurs for a very reasonable value of Q , 



1 GeV. It would 



not be difficult to write interpolation formulas, provided one tackles first the two remaining questions: the 
divergence of as{Q^) for ~ A'^, and what happens to the soft, or Pomeron contribution to cr^*(-Q2)p, 
numerically dominant when = and disappearing for ^ A^. 

In what regards the matter of the divergence of as{Q^) an interesting possibility is that, as has been 
suggested recentlyl^^', as{Q^) freezes at low Q^. Specifically, if this idea is correct, one would have an 
expression for Us like 

2 127r 
"'^•^ ^ " (33-2n/)log[(g2+M2)/yl2] ^^-^^ 
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with Af ^ a parameter that one expeets to be of the order of a typical hadronic mass (some estimates give 
M ~ 0.96 GeV). With respect to the remains of the Pomeron at large one may return to its interpretation 
as the probability of finding a real rho inside a photonl^^l and parametrize it with something like 



2ml + 



or perhaps adscribe the Pomeron to a diffractive effect and use instead an exponential form-factor: 



We will say no more about this, leaving the matter for the announced future publication. 
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5 -VERY HIGH MOMENTA AND CONNECTION WITH THE BFKL CONJECTURE 



The questions discussed in the previous section affect mostly the medium and low behaviour of 
structure functions. We will now discuss a very important matter that becomes relevant at ultra high 
energies. 

According to our equations we have (Table IV), 

F2(x',02) i?sK(g2)]-'i+(i+^)ar-\ (5.1a) 

A = 0.38 ± 0.01, rf+(l + A) = 2.406 (n/ = 4). (5.1b) 
On the other hand, the so-called BFKL equations'^^1 imply 



F^{xM') ^ a;o = ^^^, C^ = 3. (5.2) 

Q2-.00 

Eqs. (5.1) and (5.2) are, on the face of it, incompatible. What is more, while (as we have seen) (5.1) agrees 

very well with experiment, (5.2) deviates from experimental data by dozens of standard deviations. Since 
eq.(5.2) is only valid in the leading approximation and, at least to the author's knowledge, there is no control 
of the subleading corrections, one may be tempted to conclude that it should be rejected. 

The situation, however, is not as clear as may appear at first sight. First of all, the behaviour of the 
anomalous dimension matrix D(n) for large n imply that, for any fixed x the structure functions must tend 
to zero for very large Q^, contrarily to the behaviour following from eqs. (5.1) which make them grow with 
. Therefore, the behaviour (5.1) must stop at ultra large x-dependent values > Q'^{x) where F2(.t, Q^) 
should start decreasing, and perhaps turn into something like (5.2). To be precise, the conditions of vaUdity 
of (5.1) (for whose proof we refer to Sect. 6) are as follows. Assume that, at a given Qq sufficiently large for 
perturbation theory to apply, we have 

F2{x,Ql) = f{Ql)x-^''^Qo) + o{x-^^), M < Ao(Q^). 

Then, for all tQ^ > Q^, have that, for sufficiently small x 

F2{x, Q2) = (Constant)[a,(Q2)]-'^+(i+^) + ©(x""), A = Xo{Ql) independent of Q^. 

What the rcsidts of rcfs.5,6 do not say is starting from which x = x{Q'^) is this behaviour valid. The analysis 
of experimental data carried so far indicates that, when < 350 GeV^, the behaviour (5.1) holds for 
x < 10~^. In the remainder of this section we will address the question of what happens above such values 
of x and . 

To clarify the situation, and make more quantitative the analysis, we will consider the momentum sum 
rule [eq.(2.2)]. Separating the quark and gluon contributions one has (see e.g. ref.l), 



hm [\xF2{x,Q') = {el)-^^. (5.3) 



If we saturate the sum rule with the asymptotic behaviour we find that (5.2) is compatible with (5.3) provided 
C2 = (eg)[3n//(16 + 3n/)]; but if we substitute (5.1), assuming it to be valid up to a; = xq, we get 

£"dxF2{x,Q') ^ ^^^[a,(Q2)]-<^+(i+^) (5.4) 



which violates (5.3) for Q"^ > Qii^{xo) with Q^mC^o) defined by 



Bsxl^^ 16 + 3n/ 



i/d+ 
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For n/ = 4 and with the values of Bs, A we have found, say 

Bs = 10"^, A = 0.38; A = 0.10 
we get that sizeable (~ a few percent) corrections to (5.1) must occur when x ^ 10~^ already for 



^L^lOOGeV^ 



This means that the precision of the more recent HERA data may be able to reveal the corrections to (5.1), 
and even to discriminate against, or for, (5.2). 

In order to be more quantitative about this, we start by assuming that at a fixed Qq one had, exactly. 



Fs(x,Q2) = BsK(Q2)]-'i+(i+A)a;- 



(5.5) 



Then the moments at Qq are 



Ms(n,(3o) = 



Bs 



-d+(l+A) 



n-l-A' 

and, in view of the evolution equations (2.6) we get the result, valid at all Q^, 



B 



s 



n-l-X 



K(Qo)] 



-d+(l+A) 



d+{n) 



(5.6) 



an evaluation that holds to corrections of relative order 
with the momentum sum rule in that, for n = 2, (i+(2) = and thus 



In particular, (5.6) is certainly compatible 



MS (2 



Jo 



dx Fs constant. 



Next we get the behaviour implied by (5.6) for x 0. The leading behaviour is still given by eqs. 
(2.8) [or (5.1)]. Because, however now (5.6) is assumed to be exact, we may find the corrections. These are 
dominated by the first singularity of d+{n), which occurs at n = 1. After a simple calculation we find. 



Fsix,Q^) = BsiMQ^T'^^^'^^^^-^ - ^i^,Q^) + Oix), 



(5.7a) 



and A{x,Q^) is such that, for n = 1 + e, e ^ 0, it is given by the leading singularity of d+{l + e) 
[36/(33 -2n/)e]: 



[36/(33-2Ti/)e] 



= / dxx''-'^A{x,Q^). 
Jo 



At large and small x we have the asymptotic expression for the correction, 
A{x,Q^) _ 



Bs 



9iog[(ioggV.i2)/(iogg§M2)] 



(33 — 2n/)7r| loga;| 



exp . 



144| log.T| 



(33 - 2n/) 



log 



logQ2 



loggg 



(5.7b) 



Wc see how the momentum sum rule and the behaviour as 



Bs[as{Q^ 



-d+(l+A), 



get reconciled. For x ^ 0, the piece 



dominates over A{x,Q ); but, for fixed x, A dominates when oo and in fact 



(5.7) cease to be valid. It is to be noted, however, that this mechanism excludes the BFKL conjecture (5.2) 
which is never attained. 

We may generalize this last result, as the analysis depends only on the first singularity of D(n), located 
at n = 1, occurring not after the exponent of the first correction to (5.5). So we get essentially the same 
result if, at a given Q"^ one only assumes 



Fs{x, Ql) = Bs[as{Ql)]-''+^'+^^ x"^ + 0(xO). 
16 



(5.8) 



The conditions that one has the behaviour (5.2) then arc that the correction to (5.1) be of the form 
(const. )x~'^, /Lt > 0. In fact, it is not difficult to get convinced that would need an infinite set of terms 
so that one had, 

Ao = A > Ai > A2 > . . . > 0, 



A 



(5.9a) 
(5.9b) 



x=1.3xlO" 
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0.5 





10 



20 



Ve'(GeV) 



Figure 6a.- Comparison of theory and 
experiment for x = 0.013. 



x=2.0xlO 



1.0 



0.5 




10 



20 



Ve'(GeV) 



Figure 6b.- Comparison of theory and 
experiment for a; = 0.02. 



I 



Aj is given in terms of Bsj, Xj 



by a formula like (5.7b). The Bsj, Xj are parameters, not given by 
the present analysis beyond the condition Bsi > which follows from positivity of F2. Choosing them 
appropriately it would not be difficult to obtain that, although for every fixed only one term dominates, 
the envelope of (4.11) when both ^ 00, .t ^ is (5.2). 

In order to discriminate between the two possibilities, we note that, when, for a fixed x, the subleading 
terms become sizeable, they are negative, if they are like in eq. (5.7); while if they are as in (4.11) they will 
start by giving a positive contribution (imtil, in the end, the term A also here starts being important). We 
may thus employ the following strategy. We fit data points including slightly larger values of x than before. 
To be precise, we choose to add to the fit the values of ^2(3^, Q^) with a; = 1.3 x 10~^ and still we take 

12 GeV^ <Q^ < 350 GeV^. 

We keep A = 0.38, yl = 0.10 GeV fixed. We get a reasonably good fit, with xVd.o.f.= 85.1/(73 - 2), 
Bs = 2.14 X 10~^, Bjv5 = 0.54. The comparison of the result, for a; = 1.3 x 10~^ and experiment is shown 
in Fig. 5a where we have added two points with = 500 and 650 GeV^ not included in the fit; and, for 
a; = 2.0 X 10~^ in Fig. 5b, whose experimental points, however, were also not included in the fit. 

Although the results are not conclusive (one should have had more care in the treatment of other 
subleading effects, in particular those mentioned in the previous section) it would appear that data lie below 
theory, thus suggesting (5.7)* and disfavouring the BFKL conjecture, at least at the values of x, attained 
at HERA. 



* In fact one may fit the difference between the experimental points and the dominating contribution by 
a formula like (5.7b), although one should not attach too much meanig to this. 
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6. PROOF OF EQS.(2.8), (5.1), SUBLEADING CORRECTIONS AND DISCUSSION OF AL- 
TERNATIVES 



Wc repeat here the proof of the asymptotic behaviour (2.8), (5.1) for ease of reference and with a view 
to the study of subleading corrections, and to discuss the alternatives. Our assumption is that at a certain, 
fixed Qq sufficiently large for perturbation theory to be valid we have 



(6.1) 



for i = S, G. Taking moments with continuous index n as in eq.(2.5) we find that the first singularity of the 
Mi(^) Qo) functions of n occur when the integrals diverge, i.e., for n = 1 + Ai((5g). Conversely, if for any 
other we have moments diverging, 



n ■ 



A,: - 1 



then necessarily one must have 



Fi{x,Q'') ^ fi{Q^)x- 

X— >0 



(6.2) 



Because of the evolution equations for the moments. 



IJ's{n,Q^) 
l^cin, Q^) 



D(n) 



Hs{n, Ql) 



(6.3) 



the singularities of the fii{n,Q^) must originate either in those of fii{n,Q^) or in those of the anomalous 
dimension matrix D(n). By inspection, the last is seen to be analytic to the right of n = 1: so the leading 
singularity of the /Xi(n, Q"^) is identical to that of /Xi(n, Qq). Thus we may write 



Fi{x,Q^) o^^MQ^x- 

x—*0 



(6.4) 



with Aj = Xi{Qo) and thus independent of Q^. The rest of the formulas (2.8) (say) follow by diagonalizing 
the matrix D(n) and thus reducing the case to the evolution of a combination of the structure functions: 
actually, the coefficient relating Fs and Fg in (2.8f) is just the corresponding element of this diagonalizing 
matrix. 

Let us next turn to the corrections. Subtracting from e.g. Fs the leading singularity we get 



Fs{x,Q') = BsiasiQ')] 



2Nid+(l+A) ^-A 



x-^+F^s'\x'Q")- 



We can repeat the analysis for Fg''' provided it behaves, at a fixed , as ~ x~ , A(i) > 0. In this way we 
get a series like that in (5.9). If, however, the first singularity of /ii(n, Qq) (once removed that at n = 1 + A) 
lies to the left of n 1, then the dominating singularity is that of D(ri), and we get the result reported in 
(5.7). 

From the preceding analysis it follows that A, Bg cannot have a perturbative dependence on Q^; but 
they may, and likely will, have an indirect dependence via a dependence on the number of flavours excited 
at a given Q^. This has been shown to occur by Witten^-'^^l in the case of DIS on photon targets, where one 
can calculate A, and one finds a slight dependence of A on nj. This dependence (quite gcmerally now) comes 
about as follows. As discussed before, the behaviour of the structure functions is related to the singularities 
of the moments at a given Q^. Because of the operator product expansion, one has 



/x(n,^ 



fo)=Cn{p\OniQl)\p), 



(6.5a) 



and 



n{n,Q'^) = C„ 



-,D(«) 



(6.5b) 
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where Cn are the Wilson coofBcents, and the On{Qo) combinations of quark and gluon operators, renormal- 
ized at Qq, that we rather sketchily write as 

On ^ Gd^^^^^G; qf dj^^^qf , 

n n 

and / — u,d, s,c . . . runs over the varies flavours excited. By inspection, the singularities of C„, D(n) are 
seen to start at n = 1, both at leading and next to leading order. [For the singlet case; for the nonsinglet 
the relevant quantity is Dii{n) whose first singularity lies at n = 0. We send to refs. 5,6 for the details of 
study of this case]. Therefore the leading singularities are those of the operator matrix elements (p|0„|p), 
depending on the number of quark flavours that intervene. Because the more operators intervene, the more 
poles one is likely to get, this analysis also shows that one must have 

A(n/ = 5) > A(n/ 4) > A(n/ = 3) > A(n/ = 2), 

an expectation confirmed by our findings in this work. 

This argument also shows that the number of flavours to be consider should depend mostly on (and 
not on the energy, v). In fact, the integral that gives the moments in terms of the structure function recives 
contributions from all x, and not only from a; ~ 0. 

To finish this section we will discuss what happens if the assumption (6.1) fails. That is to say, if 
at all < Qq with Qq a momentum at which one may join the nonperturbative Pomeron regime, and 
perturbative QCD (assumed to exist) one has 

F,{x,Ql)^0{x''). 

This was discussed long ago by De Riijula et al.^° and has been revived recently (see e.g. ref. 21). As 

explained above, this means that the singularities in n of the (p|On(Qo)!p) lie at, or to the left of n = 1 
[else we would have had (6.1)]. Then the behaviour of the Fj will be dictated by the singularities of D(n), 
for i = S, G or Z)ii(n) for the nonsinglet case. The analysis is essentially like for the corrections A to F2 
discussed in the previous section and we immediately obtain, 

F2{x, g2) ^ exp y^6(log(logg2/logQ2))|ioga;|, (6.6a) 

and, for the nonsinglet, 

FNsix, Q^) ~ exp ^6'(log(loggVlogQ2))|iogx|. (6.6b) 

The constants b, b', and the more detailed form of the behaviour depend on the assumption one makes at 
= Qoi see e.g. ref. 21. For example, if we assume that Fs{x,Qq) ~ C, then the formula for Fs{x,Qq) 
becomes like (5.7b), with the obvious changes and in particular, 

b- "^^^ 
~ 33 - 2n/" 

It is perhaps not idle to emphasize that (6.6) and (6.1) are mutually exclusive. If, for any one has 

(6.1) with A strictly positive, one necessarily gets (2.8) and, conversely. Only experiment may discriminate 
between the two alternatives since the value of A is not fixed by perturbative QCD. Thus either (2.8) or 
(6.6), but not both, will fit experiment if the last is sufiiciently precise. 

The fits with (2.8) have been described in the previous sections. In what regards eqs.(6.6), I have found 
it impossible to fit actual experimental data with them. Eyeball inspection shows that alkhe NS structure 
functions behave as x^ near a; = 0; and evey nonbiased fit, including those in refs. 10, 11 (not to mention 
ours in the present paper!) give an F2 behaving also as a power x~^, with similar exponents A for all 
(between 0.23 and 0.40) in the region of HERA data. On the other hand, any fit with eqs. (6.6), performed 
fixing x to get rid of the large uncertainties in (6.6), give impossibly large values of F2 as grows beyond 
100 GeV^ or impossibly smafl for < 40GeV^ 
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